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ABSTRACT 

Seismic  Data  Collection  System  (SDCS)  stations  were  deployed  at  the 
Nevada  Test  Site  (NTS)  and  two  in  the  Eastern  United  States  (EUS)  to  measure 
magnitude  residuals  and  spectral  differences  between  NTS  and  EUS  stations.  The 
deployment  was  intended  to  determine  the  degree  of  anelastic  attenuation  under 
NTS.  At  the  Climax  Stock  (0B2NV)  station,  the  teleseismic  body-wave  magnitudes 
are  .17  magnitude  units  (m.u.)  lower  than  at  the  EUS  stations.  The  magnitudes 
at  the  two  Pahute  Mesa  sites  are  about  .2  m.u.  higher  than  at  0B2NV,  a  difference 
that  can  be  accounted  for  by  amplification  of  low  velocity  volcanics  under  Pahute 
Mesa.  Thus  local  geology  must  be  taken  into  account  in  order  to  evaluate  the 
attenuation  under  a  given  site  using  amplitude  data.  This  correction  shows  that 
with  respect  to  corrected  amplitude  levels  0B2NV  is  equivalent  to  Pahute  Mesa 
stations.  At  all  NTS  sites  the  higher  frequency  content  of  P  waves  is  signi¬ 
ficantly  less  than  in  the  EUS.  Available  data  suggests  an  attenuation  effect 
of  about  .2  m.u.  under  NTS.  Measurements  at  a  few  other  Western  United  States 
(WUS)  sites,  including  the  site  of  the  SHOAL  explosion,  yield  similar  figures. 
These  results  confirm  that  body-waves  suffer  considerable  anelastic  losses  traver¬ 
sing  the  mantle  under  the  WUS,  including  the  NTS  sites. 
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INTRODUCTION 


Many  studies  have  confirmed  the  existence  of  the  Western  United  States 
(WUS)  magnitude  anomaly  for  short-period  teleseismic  P  arrivals  that  was  dis¬ 
covered  in  the  early  sixties  (Guyton,  1964;  Evernden  and  Clark,  1970;  Booth, 
Marshall  and  Young,  1974;  Der,  Masse  and  Curski,  1975;  North,  1976),  The 
anomaly  demonstrates  its  presence  in  the  form  of  a  regional  reduction  of  tele- 
seismic  magnitude  measurements  at  WUS  receiving  stations  compared  to  those 
located  in  the  Eastern  United  States  (EUS).  A  similar  anomaly  seems  to  exist 
for  short-period  S  waves  (Der,  Masse  and  Curski,  1975).  In  addition,  spectral 
differences  exist  between  short-period  body-waves  crossing  the  upper  mantle 
under  the  WUS  and  those  crossing  the  mantle  under  shields  or  stable  platforms 
(Der  and  McElfresh,  1977). 

Most  seismologists  attribute  the  regional  magnitude  anomaly  and  spectral 
differences  to  anelastic  attenuation  in  the  mantle  under  the  WUS.  Other  tvpes 
of  geophysical  anomalies  in  the  WUS,  which  indicate  partial  melting  in  the 
mantle  under  the  WUS,  tend  to  support  this  view.  This  partial  melting  of 


Guyton,  J.  W.  (1964).  Systematic  deviations  of  magnitude  from  hody-waves  at 
seismograph  stations  in  the  United  States.  Proceedings  of  VESIAC  Con¬ 
ference  on  Seismic  Event  Magnitude  Determination,  University  of  Michigan, 
441 0-71-X. 

Evernden,  J.  and  D.  M.  Clark  (1970).  Study  of  teleseismic  P.  TI  Amplitude 
data,  Phys.  Earth.  Plant  Int.  ,  4_,  24-31. 

Booth,  D.  C. ,  P.  D.  Marshall,  and  J.  B.  Young  (1974).  Long-  and  short-period 
amplitudes  from  earthquakes  in  the  range  0°-114°,  Geophvs.  J.  R.  Astr. 
Soc. ,  39,  523-538. 

Der,  Z.  A.,  R.  P.  Masse  and  J.  P.  Curski  (1975).  Regional  attenuation  of 
short-period  P  and  S  waves  in  the  United  States,  Geophys.  J.  R.  Astr. 

Soc. ,  40,  85-106. 

North,  R.  G.  (1976).  Station  biases  in  body-wave  magnitude  (abstract). 
Transaction,  American  Geophysical  Union  (EOS) ,  57 ,  955. 

Der,  7.  A.  and  T.  W.  McElfresh  (1977).  The  relationship  between  anelastic 

attenuation  and  regional  amplitude  anomalies  of  short-period  P  waves  in 
North  America,  Bull.  Seism.  Soc.  Am.,  67,  No.  5,  1303-1317. 
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material,  according  to  experimental  and  theoretical  studies  (Walsh,  1968,  1969), 
must  lead  to  the  increase  of  anelastic  attenuation.  The  spectral  differences 
of  P  wave  arrivals  in  the  WUS  versus  those  of  the  EUS  are  such  that  the  size 
of  the  magnitude  anomaly  is  roughly  consistent  with  spectral  differences, 
assuming  a  simple  exponential  attenuation  factor  of  the  form  exp (-irf At*)  where 
At*  is  derived  from  slopes  of  spectral  ratios. 

Evoking  reciprocity,  there  is  the  implication  that  anelastic  attenuation 
causes  a  bias  in  the  empirical  magnitude-yield  relationships  based  upon  NTS 
data.  Therefore,  determining  the  magnitude  of  the  attenuation  effect  is  im¬ 
portant. 

Previous  studies  have  established  the  size  of  the  average  EUS-WUS  regional 
magnitude  bias  and  spectral  differences  fairly  well  (Am^  ~  .25,  At*  -  .2). 

This  cannot  be  said  however,  for  individual  receiving  stations,  because  both 
the  Am^  and  At*  can  differ  from  the  regional  value  at  individual  sites.  To 
estimate  the  amount  of  anelastic  attenuation  under  individual  stations,  a 
detailed  study  of  both  Amfa  and  At*  must  be  made  using  many  teleseismic  events. 
Because  crustal  structures  can  also  affect  both  quantities,  efforts  must  be 
made  to  estimate  and  correct  for  the  crustal  effects  using  the  best  estimate 
of  the  local  crustal  structure.  Such  corrections  for  crustal  structure  are 
difficult  to  make  because  the  structure  is  not  sufficiently  known  at  most  places. 
The  properties  of  near  surface  layers,  which  have  the  most  effect  on  both  Ara^ 
and  At*,  are  especially  difficult  to  estimate.  Thus,  regional  geological  infor¬ 
mation  is  often  not  sufficient  and  local  surveys  may  be  necessary. 

The  authors'  past  work  in  simulating  the  effects  of  crustal  structure  in¬ 
dicated  that  Am^  may  be  influenced  more  by  crustal  structure  than  by  typical 
differences  in  At*.  On  the  other  hand  formulating  a  structure  that  would  be 
able  to  remove  all  high  frequencies  from  the  spectrum  as  anelastic  attenuation 
does  is  difficult  if  not  impossible. 


Walsh,  J.  B.  (1968).  Attenuation  in  partially  melted  material,  J.  Geophys. 
Res. ,  73,  2209-2216. 

Walsh,  J.  B.  (1969).  New  analysis  of  attenuation  in  partially  melted  rock, 
J.  Geophys.  Res.,  74,  4333-4337. 
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To  investigate  the  question  of  the  existence  of  high  anelastic  attenua¬ 
tion  under  NTS,  three  Special  Data  Collection  System  (SDCS)  seismic  stations 
were  deployed  to  that  area.  One  station,  0B2NV,  was  placed  on  Climax  Stock, 
a  Cretaceous  granitic  intrusive  body,  the  site  of  the  HARDHAT  and  PILE-DRIVER 
nuclear  explosions.  Two  other  stations,  NTNV  and  NT2NV,  were  placed  on  the 
Pahute  Mesa,  a  deep  volcanic  caldera  (Silent  Canyon  Volcanic  Crater) ,  which 
is  filled  with  low  velocity  [a  -  3.6  km/sec  average]  Tertiary  volcanic  tuffs 
and  rhyolites  to  a  depth  of  at  least  4  km.  Spence  (1974)  predicted,  based 
upon  gravity  data,  a  thickness  of  about  5  km.  Two  other  SDCS  stations  remained 
at  their  previous  locations:  Red  Lake,  Ontario  (RKON)  and  Houlton,  Maine 
(HNME) .  A  key  station  at  NTS  is  0B2NV  which,  like  RKON,  is  located  on  granite, 
the  test  medium  of  greatest  interest.  HNME  is  located  in  a  region  where  there 

It 

are  indications  of  some  attenuation  in  the  upper  mantle  (Solomon  and  Toksoz, 

* 

1970;  Der,  Masse  and  Gurski,  1975).  NT2NV  and  NTNV  are  less  diagnostic  because 
they  are  located  above  the  thick  low-velocity  caldera  fill  and  (possibly)  a 
deep,  buried  volcanic  plug  which  was  hypothesized  by  Spence  (1974).  The  loca¬ 
tions  of  the  stations  used  in  this  study  are  shown  in  Figure  1.  In  this  figure 
we  also  show  the  LRSM  stations  SZNV,  SEMN,  CUNZ,  and  EKNV  also  analyzed  in  this 
report. 


Spence,  W.  (1974).  P  wave  residual  differences  and  inferences  on  an  upper 
mantle  source  for  the  Silent  Canyon  Volcanic  Centre,  Southern  Great 
Basin,  Nevada,  Geophys.  J.  R.  Astr.  Soc.,  38,  505-523. 

Solomon,  S.  C.  and  M.  N.  Toksoz  (1970).  Lateral  variation  of  attenuation  of 
P  and  S  waves  beneath  the  United  States,  Bull.  Seism.  Soc.  Am.,  60, 
819-838. 
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Figure  1.  Map  of  station  locations  used  in  this  study 


DATA  ANALYSIS 


k 


Body-Wave  Magnitudes 


Short-period  data  at  the  SDCS  sites  were  recorded  on  digital  tape  at  20 
samples/second.  Amplitudes,  dominant  periods  and  arrival  times  were  read  for 
each  available  station  using  the  CRT  display  of  traces  on  the  PDP-15  computer 
for  a  set  of  119  events  that  covered  the  period  between  1  September  1976  to 
31  March  1977. 


Events  used,  along  with  the  amplitude  and  period  readings  and  magnitudes 
computed  at  each  station,  are  given  in  Appendix  A.  Magnitudes  are  also  given 
in  the  table  for  each  event,  when  available.  Magnitudes  (m^)  were  computed  with 
the  formula 

"H)  =  logio(P  +  B(A)  (1) 

where  A  is  the  amplitude  in  millimicrons  (my)  (corrected  for  instrument  ampli¬ 
fication  at  the  measured  period)  and  T  is  the  dominant  period  of  the  signal. 

B(A)  are  Gutenberg  and  Richter's  (1956)  (peak-to-peak)  distance  correction 
factors  for  aero  hypocentral  depth.  We  also  computed  "magnitudes"  without 
dividing  by  the  period,  but  correcting  the  amplitudes  for  instrument  response 
at  the  measured  period.  We  shall  call  these  values  "uncorrected"  in  this  report. 

(uncorrected)  =  )  +  B(A)  (2) 

b  1U  my 

The  "uncorrected"  magnitudes  minimize  the  magnitude  effect  due  to  attenu¬ 
ation.  As  the  signal  decreases  in  amplitude  due  to  attenuation,  the  lower  fre¬ 
quencies  are  less  affected  and  the  dominant  period  increases.  The  system  mag¬ 
nification  is  less  at  long  periods  so  that  the  amplitude  corrected  for  instru¬ 
ment  response  at  the  larger  value  of  T  tends  to  remain  constant.  Dividing  by 
T  tends  to  again  reduce  the  magnitude.  Thus  when  T  is  not  divided  out,  the 
magnitude  tends  to  remain  the  same  as  attenuation  increases.  The  greatest 
difference  in  amplitude  is  seen  on  the  raw  traces,  uncorrected  for  instru¬ 
ment  response.  These  differences  are  not  given  in  this  report. 

Results  of  this  study  are  displayed  in  a  set  of  histograms  showing 
magnitude  differentials,  dominant  periods  and  At*  between  all  possible  pairs 


Gutenberg,  B.  and  C.  F.  Richter  1956,  Magnitude  and  energy  of  earthquakes, 
Ann.  Geof.  (Rome),  1-15. 


-13- 


T" 


of  individual  stations.  In  each  histogram  we  calculated  the  standard 

deviation  (a),  the  mean  value,  and  its  approximate  95%  confidence  limits 

±2a  where  a  =  o//n  (assuming  normality).  For  small  samples  (n  <  30) 
mean  mean 

we  give  confidence  limits  based  on  the  "Student's  t"  distribution  (Hoel, 

1954).  We  also  print  the  total  number  of  data  points  used  on  each  histo¬ 
gram.  We  did  not  use  any  readings  at  epicentral  distances  greater  than  85°. 
This  was  done  to  eliminate  the  effect  of  apparent  low  Q  due  to  diffraction 
close  to  the  core-mantle  boundary. 

Figures  2  through  29  show  histograms  for  magnitude  differentials,  both 
corrected  and  uncorrected,  as  well  as  differences  in  dominant  P  wave  periods 
for  all  possible  station  pairs.  Period  histograms  for  HNME  were  not  included 
because  a  different  instrument  response  at  this  station  made  any  comparison 
invalid. 

Figure  2,  which  shows  the  RK0N-0B2NV  pair,  is  the  most  diagnostic  dis¬ 
cussed  here.  The  corrected  magnitudes  are  significantly  less  at  0B2NV  than 
at  RKON.  Despite  the  large  number  of  data  points  (76)  used,  the  confidence 
limits  are  still  quite  wide  (±.114).  Although  the  sign  of  the  mean  is  still 
positive,  which  may  indicate  lower  amplitudes  at  0B2NV,  the  "uncorrected" 
magnitude  differentials  (Figure  3)  do  not  show  a  significant  bias.  While 
this  situation  does  not  rule  out  the  existence  of  an  uncorrected  magnitude 
(amplitude)  differential,  it  does  show  that  the  differential  can  be  anywhere 
between  -.04  to  .21  (95%  limits).  Since  about  seven  months  of  data  are  shown, 
and  the  scatter  yields  95%  confidence  limits  of  +.1,  the  sites  should  be 
occupied  for  about  two  years  with  the  same  average  rate  of  data  acquisition, 
quadrupling  the  data  to  narrow  the  confidence  limits  to  about  ±.05  magnitude 
units.  The  difference  in  average  P  wave  periods  (Figure  4)  is  also  signifi¬ 
cant  which,  since  the  dominant  periods  at  this  station  are  longer,  supports 
the  hypothesis  of  higher  attenuation  under  0B2NV.  (Average  P  wave  period  at 
RKON  is  .70  seconds  vs  0.88.) 

In  contrast,  when  NT2NV  and  NTNV  are  compared,  no  significant  difference 
is  found  between  the  two  stations  either  in  magnitude  or  dominant  period 


Hoel,  P.  0.  1954,  Introduction  to  Mathematical  Statistics:  New  York, 
John  Wiley  &  Sons. 
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Figure 


3.  Histogram  of  uncorrected  magnitude  differentials  RK0N-0B2NV 
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4.  Histogram  of  differences  in  dominant  P  wave  periods  RKON-OB2NV 


CORRECTED  A  mL 

D 

Figure  5.  Histogram  of  magnitude  differentials  NTNV-NT2NV 
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Figure  8.  Histogram  of  magnitude  differentials  r!TNV-0B2*JV 
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Figure  9. 


Histogram  of  uncorrected  magnitude  differentials  VTll\’-OP2W 
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12.  Histogram  of  uncorrected  magnitude  differentials  NT2NV-OB2NV 
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Figure  13. 


Histogram  of  differences  in  dominant  P  wave  periods  NT2NV-OB2NV 
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Figure  14.  Histogram  of  magnitude  differentials  RKON-NTNV 


NUMBER  NUMBER 


Figure  17. 


Flistogram  of  magnitude  differentials  RK0N-NT2NV 
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'igure  70.  Histogram  of  magnitude  differentials  HNME-RKON 


Figure  21. 


Histogram  of  uncorrected  magnitude  differentials  HKMg-RKOlJ 
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rigure  24.  Histogram  of  magnitude  differentials  FNME-NT2NV 
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Figure  25.  Histogram  of  uncorrected  magnitude  differentials  H!CME-NT2NV 
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Figure  26.  Histogram  of  magnitude  differentials  HNME-0R2NV 


Figure  27  Histogram  of  uncorrected  magnitude  differentials  HNME-OB2NV. 
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Figure  28.  Histogram  of  List-RKON  magnitude  differentials 
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Figure  29.  Histogram  of  List-0B2NV  magnitude  differentials 
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(Figures  5  through  7). 

Considerable  difference  exists,  however,  between  0B2NV  and  the  two  Pahute 
Mesa  stations  (Figures  8  through  18).  Amplitudes  are  considerably  higher  and 
the  dominant  periods  shorter  at  Pahute  Mesa  than  on  the  Climax  Stock. 

Magnitudes,  and  amplitudes,  at  Pahute  and  RKON  are  roughly  on  the  same 
level  and  the  dominant  periods  are  intermediate  between  RKON  and  0B2NV  (Fig¬ 
ures  14  through  19).  Calculations  in  the  subsequent  sections  of  this  report 
shows  that  because  of  the  low— velocity  crustal  layers,  amplitudes  at  1  anute 
are  0.2  log  amplitude  higher  than  they  would  be  on  granite,  a  result  that  brings 
the  RK0N-NT2NV  and  RKON-NTNV  amplitude  differences  into  rough  agreement  with 
the  RK0N-0B2NV  differences. 

Note  that  the  mean  of  standard  deviations  of  Am,  measurements  at  the 

r> 

three  NTS  sites,  compared  to  RKON,  is  almost  twice  the  mean  standard  deviation 
of  the  relative  Am^  among  pairs  of  NTS  stations.  This  situation  is  at  least 
partially  a  result  of  the  events'  radiation  patterns,  which  cause  a  greater 
variability  between  azimuths  to  the  RKON-NTS  station  pairs  from  the  epicenters, 
while  the  NTS  stations  are  set  at  almost  the  same  point  in  the  radiation  pat¬ 
terns,  suggesting  an  alternative  way  to  determine  RKON  to  NTS  bias  using  a 
reference  station,  denoted  by  X,  close  to  NTS  that,  together  with  RKON,  has 
been  operational  for  a  long  time.  Thus,  the  bias  of  station  X,  relative  to 
RKON,  could,  using  a  large  number  of  events,  be  measured  with  great  precision. 
From  this  point  deducing  the  station  X-NTS  bias  (and  RKON-NTS)  using  fewer 
events  since  this  pair  would  be  less  affected  by  radiation  patterns  would  be 
a  simple  matter.  The  proximity  of  station  X  to  NTS  could  also  be  useful  in 
measuring  identical  portions  of  the  wavetrain,  which  has  been  done  for  the 
SZNV-MNNV  station  pair  in  the  past  (Der  and  McElfresh,  1975).  The  SRO  sta¬ 
tion  ANMO  may  also  be  used  this  way  because  North  (1976)  determined  its 
magnitude  bias  with  great  precision. 

Der,  Z.  A.  and  T.  W.  McFlfresh  (1975).  Short-period  P  wave  attenuation  along 
various  paths  in  North  America  as  determined  from  P  wave  spectra  of  the 
SAI.MON  nuclear  explosion,  SDAC-TR-75-16 ,  Teledyne  Geotech,  Alexandria, 
Virginia  22314. 

North,  R.  G.  (1976).  Station  biases  in  bodv-wave  magnitude  (abstract). 
Transactions,  American  Geophysical  Fnion  (E$S),  57,  955. 
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Located  in  an  area  of  high  microseismic  noise  and  suspected  higher 
attenuation,  HNME,  the  last  station  discussed  here,  is  the  least  valuable  for 
our  study  (Solomon  and  Toksoz,  1970;  Der,  Masse  and  Gurski,  1975).  Because 
of  the  high  noise  level  there  were  considerably  fewer  readings  at  this  station, 
but  no  significant  bias  was  found  relative  to  RKON  (Figures  20  and  21).  Com¬ 
parisons  of  dominant  periods  involving  HNME  readings  are  not  shown  because 
instrument  response  at  this  station  was  significantly  less  peaked  at  high 
frequencies  than  at  other  stations  used  in  this  study.  The  magnitude  compari¬ 
sons  are  valid,  however,  since  instrument  response  is  corrected.  Comparison 
of  HNME  magnitudes  with  the  NTS  stations  are  inconclusive  due  to  wide  confi¬ 
dence  limits  (Figures  22  through  27). 

Comparison  of  magnitudes  in  the  list  with  our  data  show  that  those 
magnitudes  agree  with  readings  at  0R2NV  and  are,  consequently,  significantly 
lower  than  RKON  magnitudes  (Figures  28  and  29). 

Table  I  summarizes  the  results  of  the  magnitude  study  at  SDCS  stations. 
Spectral  Ratio  Studies 

In  addition  to  magnitude  measurements  relative  t*  values  were  also 
computed  for  all  combinations  of  station  pairs  for  a  set  of  selected  events. 

The  At*  estimates  were  obtained  by  fitting  a  straight  line  to  the  logarithm 
of  the  amplitude  spectral  ratio  between  each  pair  of  stations.  The  signal 
window  was  of  12  sec  duration.  This  procedure  is  identical  to  that  described 
in  a  previous  report  (Der  and  McElfresh,  1976b),  but  for  completeness  it  is 
summarized  in  the  next  section.  A  3-second  signal  window  was  also  used  for 
many  of  the  events  because  we  suspected  that  the  early  part  of  the  signal 
might  be  more  diagnostic.  However,  the  basic  features  of  most  spectra  com¬ 
puted  from  3-second  signal  windows  remained  the  same  as  those  from  12-second 
windows,  although  the  variance  of  At*  increased.  (In  this  report  only  the 
At*  computed  from  12-second  windows  are  used.)  The  individual  At*  determina¬ 
tions  are  tabulated  in  Table  II. 


Der,  Z.  A.  and  T.  W.  McElfresh  (1976b).  The  effect  of  attenuation  on  the 
spectra  of  P  waves  from  nuclear  explosions  in  North  America,  SDAC-TR- 
76-7,  Teledyne  Geotech,  Alexandria,  Virginia  22314. 


TABLE  I 
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Summary  of  SDCS  Magnitude  Study 


STATION  PAIR  CORRECTED  UNCORRECTED  mb  PERIODS 


- 

U 

2ct 

u 

u 

2a 

U 

u 

2a 

u 

RKON-OB2NV 

0.184 

0.114 

0.085 

0.123 

-0.185 

0.048 

-  - 

4 

h 

NTNV-NT2NV 

0.017 

0.068 

0.017 

0.077 

0.  000 

0.054 

> 

NTNV-OB2NV 

0.240 

0.099 

0.182 

0.107 

-0.099 

0.061 

V 

NT2NV-OB2NV 

0.238 

0.076 

0.177 

0.089 

-0.109 

0.052 

1 

RKON-NTNV 

-0.107 

0.138 

-0.136 

0.144 

-0.065 

0.070 

- 

RKON-NT2NV 

-0.080 

0.110 

-0.113 

0.119 

-0.068 

0.052 

* 

HNME-RKON 

-0.051 

0.109 

0.030 

0.127 

f 

HN21E-NTNV 

-0.197 

0.177 

-0.130 

0.206 

* 

; 

HNME-NT2NV 

-0.077 

0.141 

-0.062 

0.161 

HNME-OB2NV 

0.138 

0.163 

0.101 

0.177 

LST-RKON 

-0.253 

0.098 

s 

i 

LST-OB2NV 

-0.076 

0.115 
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Results  of  all  At*  measurements  between  SDCS  station  pairs 
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Figures  B1  to  B20  in  the  Appendix  display  the  signal  traces  for  events 
studied.  An  "X"  marks  traces  omitted  from  the  compilation  of  spectral  results 
because  of  poor  S/N  ratio.  An  asterisk  (*)  marks  traces  excluded  because  the 
hypo-central  distance  exceeded  85°.  Figures  Cl  to  C80  show  all  spectra  com¬ 
puted  and  Figures  D1  to  D80  show  all  amplitude  spectral  ratios  with  straight- 
line  fits. 

The  resulting  mean  At*  values,  as  well  as  the  standard  deviation,  standard 
deviation  of  the  mean,  and  the  95%  confidence  limits  on  the  mean  are  displayed 
as  histograms.  Positive  values  of  At*  indicate  that  the  station  spectrum  in 
the  denominator  is  richer  in  high  frequencies.  The  general  trend  emerging 
from  these  histograms  is  similar  to  the  trend  derived  from  magnitude  and 
period  data.  0B2NV  is  again  distinguished  by  the  low-frequency  character  of 
its  P  wave  arrival  when  compared  to  RKON  (Figure  30).  The  mean  At*  is  about 
0.196,  which  can  explain  a  0.23  magnitude  bias  at  1  Hz.  The  difference 
between  NT2NV  and  NTNV  is  small  and  not  significant  (Figure  31),  but  both 
stations  are  richer  in  high  frequencies  than  0B2NV  (Figures  32  and  33).  Both 
Pahute  Mesa  stations  are,  in  turn,  significantly  lower  in  high  frequency  con¬ 
tent  of  P  waves  than  RKON  (Figures  34  and  35),  indicating  their  intermediate 
position  between  0B2NV  and  RKON.  No  conclusions  can  be  made  at  the  95%  confi¬ 
dence  level  from  the  histograms  involving  HNME  because  of  the  small  number  of 
available  ratios  (Figures  36  through  39). 

In  summary,  the  spectral  analyses  provide  a  picture  that  is  qualita¬ 
tively  consistent  with  magnitude  biases  and  visual  observations  of  dominant 
periods.  The  most  important  results  of  the  At*  study  are  shown  in  Table  III. 

Note  that  the  observed  differences  in  dominant  periods  of  P  waves  in  the 
EUS  and  WUS,  respectively,  are  of  the  order  of  .1-.2  sec  (the  average  period 
for  teleseismic  arrivals  is  about  1  sec).  This  difference  is  much  less  than 
differences  between  dominant  periods  of  NTS  explosions  and  Kazakh  explosions 
of  similar  body-wave  magnitude  observed  at  common  stations,  which  were  on  the 
average  about  .4  sec  (Springer,  personal  communication).  Simulations  performed 
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Figure  31. 


Histogram  of  At*  derived  from  the  NT2NV/NTNV  spectral  ratio 
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Figure  33. 
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Histogram  of  At*  derived  from  the  NT2NV/RK0N  spectral  i 
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Figure  35. 


Histogram  of  At*  derived  from  the  NTNV/RKON  spectral  ra 


Figure  38. 


Histogram  of  At*  derived  from  the  NTNV/HNMF  spectral  ratio 
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Figure  39. 


Histogram  of  At*  derived  from  the  0B2NV/HN>fF  spectral  ratio 


TABLE  III 


Summary  of  At*  Measurements 


STATION  PAIR 

MEAN 

95%  CONFIDENCE 
LIMITS  ON  MEAN 

1. 

NT2NV/HNME 

-0.057 

0.130 

2. 

NT2NV/RK0N 

0.091 

0.065 

3. 

NT2NV/OB2NV 

-0.059 

0.050 

4. 

NT2NV/NTNV 

0.011 

0.050 

5. 

NTNV/HNME 

0.027 

0.168 

6. 

NTNV/RKON 

0.148 

0.098 

7. 

NTNV/OB2NV 

-0.079 

0.055 

8. 

OB2NV/HNME 

0.115 

0.134 

9. 

OB2NV/RKON 

0.196 

0.055 

10. 

RKON/HNME 

-0.103 

0.117 

-39- 


on  synthetic  waveforms  (Blandford  and  Der,  1977;  Blandford,  1976)  showed  that 
the  differential  At*  'v  .2,  which  was  postulated  between  EPS  and  WPS,  could 
not  explain  a  .4  sec  difference  in  the  dominant  period,  and  the  differential 
of  dominant  periods  should  have  been  in  the  range  .1-.2  sec  for  typical  tele- 
seismic  arrivals  with  corner  frequencies  of  around  1.-1.5  Hz.  Therefore, 
most  of  the  .4  sec  differential  that  Springer  observed  must  have  been  due  to 
causes  other  than  attenuation  such  as  the  difference  in  scale  depth  of  explo¬ 
sions  at  NTS  versus  Kazakh  or  some  of  the  effects  of  shot  media.  Thus  it  is 
incorrect  to  attribute  the  EUS-WUS  magnitude  differential  to  the  dominant  per¬ 
iod  alone,  through  the  division  in  A/T.  The  EUS-WUS  magnitude  differential, 
therefore,  was  mostly  the  result  of  amplitude  differences,  which  we- e  signi- 
ficantly  lower  in  the  WUS  on  the  95%  confidence  level  (Der,  Masse  and  Gurski, 
1975).  The  contribution  of  dominant  period  to  the  differential  must  be  less 
than  or  roughly  equal  to  .1  magnitude  unit. 

Procedures  Followed  in  Computing  t* 

The  short-period  vertical  P  wave  seismograms,  including  a  noise  sample 
prior  to  the  arrival  of  the  P  waves,  were  anti-alias  filtered  and  digitized 
from  magnetic  tape  at  20  samples/sec. 

A  time  window  comprising  9  sec  of  the  signal  and  4  sec  of  the  preceding 
noise  was  tapered  with  a  Parzen  window  and  Fourier  transformed,  and  then  the 
power  spectrum  was  computed  by  multiplying  the  Fourier  spectrum  by  its  conju¬ 
gate.  This  power  spectrum  was  then  smoothed  by  a  12-point  running  average. 

The  shift  of  the  time  window  was  designed  to  avoid  heavy  tapering  of  the  first 
arrival.  A  noise  power  spectrum  ahead  of  this  window  was  computed  using  an 
identical  procedure  and  subsequently  subtracted  from  the  spectrum  of  the 
window  containing  the  signal. 


Blandford,  R.  R.  and  Z.  A.  Der  (1977).  Gross  characteristics  of  short-period 
seismic  P  waveforms  from  nuclear  explosions  as  functions  of  source  type, 
medium,  instruments,  attenuation,  and  burial  depth.  SDAC  technical 
report  in  preparation.  Teledyne  Geotech,  Alexandria,  Virginia  22314. 

Blandford,  R.  R.  (1976).  Experimental  determination  of  scaling  laws  for  con¬ 
tained  and  cratering  explosions,  SDAC,  TR-76-3,  Teledvne  Geotech, 

Alexandria,  Virginia  22314. 
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The  source  spectra,  determined  from  scaling  or  close-in  measurements, 
were  then  modified  with  the  instrument  response  and  smoothed  with  a  12-point 
running  average.  This  procedure  was  intended  to  simulate  an  identical  treat¬ 
ment  given  to  the  data  spectra  (the  effect  of  the  Parzen  window  is  not  very 
significant  relative  to  the  12-point  smoothing).  Instead  of  correcting  the 
spectra  for  pP  or  spall,  which  for  any  sizeable  surface  reflection  coefficient 
would  introduce  spurious  peaks  into  the  spectrum,  smoothing  was  relied  upon 
to  eliminate  any  modulation  present  due  to  pP  or  other  causes.  The  fit  of 
straight  lines  to  spectral  ratios  constitutes,  in  effect,  a  final  smoothing. 

We  used  only  the  portion  of  spectrum  between  .5  and  4.  Hz,  and  in  fitting 
straight  lines  to  the  spectral  ratios  the  signal  power  spectrum  was  required 
to  exceed  the  noise  by  a  factor  of  three.  Portions  of  spectra  not  satisfying 
this  criterion  were  not  used  in  the  least  squares  fit.  Furthermore,  to  avoid 
problems  with  the  dynamic  range  of  the  system,  the  portions  of  power  spectra 
which  were  down  more  than  2.5  orders  of  magnitude  from  the  peak  were  also 
disregarded. 

The  amplitude  spectral  ratio  was  then  obtained  by  taking  the  square  root 
of  the  power  ratio.  Slopes  (a)  were  measured  from  the  graphs,  and  the  values 
of  t*  were  computed  from  these  by  the  application  of  the  following  formula: 


ESTIMATION  OF  THE  EFFECT  OF  LOCAL  GEOLOGY 


Since  the  geological  structure  under  the  three  NTS  sites  is  quite 
different,  the  effects  of  geology  on  the  observed  amplitudes  must  be  consid¬ 
ered.  Stations  NT2NV  and  NTNV  are  at  Pahute  Mesa,  which  is  located  over  a 
deep  volcanic  caldera  (Silent  Canyon  caldera)  filled  with  relatively  low 
velocity  (a  =  3.6  km/sec)  tertiary  volcanics.  The  total  thickness  of  these 
volcanics,  tuffs  and  rhyolite  flows  is  not  known,  but  deep  drilling  on  the 
Pahute  Mesa  penetrated  the  caldera  fill  to  a  depth  of  about  4  km  without 
reaching  the  bottom  of  the  tertiary  volcanic  sequence.  The  0P>2NV  site,  on 
the  other  hand,  is  located  on  a  granite  stock  which  intrudes  into  high- 
velocity  paleozoic  rocks.  NT2NV  is  close  to  the  site  of  the  POOL  nuclear 
explosion  while  NTNV  is  located  near  the  KASSERI  test  site.  Some  geological 
information  is  available  from  drill  holes  at  both  nuclear  test  sites  and 
crustal  models  for  each  SDCS  station  location  were  constructed  using  this 
information.  Table  IV  gives  the  crustal  models  for  each  station. 

The  crustal  amplification  was  estimated  by  computing  the  vertical  dis¬ 
placement  for  a  P  wave  pulse  incident  on  the  models  at  three  angles  of 
incidence,  =  20,  25,  and  30  degrees  measured  from  the  vertical.  The  pulse 
used  was  obtained  from  the  von  Seggern-Blandford  source  model  with  cube-root 
scaling,  and  it  has  the  shape  of  a  far-field  displacement  from  a  50-kt 
nuclear  explosion  in  granite.  The  spectrum  of  this  pulse  resembles  a  typical 
teleseismic  arrival  with  most  of  the  energy  below  1  Hz.  Figures  40  through 
42  show  the  frequency  responses,  time  domain  waveforms,  and  readings  at  the 
three  sites. 

Table  V  displays  the  results  of  computations.  The  caldera  filling  seems 
to  account  for  about  a  .2  magnitude  increase  relative  to  0B2NV  at  Pahute  Mesa. 
Correction  for  this  effect  brings  the  amplitudes  at  Pahute  Mesa  down  to  the 
same  level  as  0B2NV.  Because  at  low  frequencies  this  amplitude  variation  is 
expected  to  vanish,  the  observed  spectrum  at  Pahute  Mesa  must  be  enhanced  at 
high  frequencies. 

In  addition  to  trying  to  assess  the  effects  of  local  geology  attempts 
were  also  made  to  find  a  pattern  in  the  azimuth-distance  distribution  magni¬ 
tude  residuals,  but  no  consistent  behavior  emerged.  Figures  43a  through  43e 
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Crustal  Models 
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3.60 

2.80 

10. 

6.  70 
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3.60 

2.80 

10. 

6.70 

4.00 

2.95 

CO 

7.80 

4.50 

3.30 

i 

1 


RECEIVER  CRUST 


NT2NV 


Figure  41.  Crustal  response  calculations  at  NT2NV 
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TABLE  V 


Summary  of  the  Differences  in  Crustal  Response 
Expressed  in  Magnitude  Units 


PULSE  AMPLITUDES  AMPLITUDE  RATIOS 

(Peak  to  Peak) 


ANGLE  OF 
INCIDENCE 

X 

J 

f - 

0B2NV 

NT2NV 

N"NV^ 

' NT2/0B 

_ , 

NT/OB 

ro 

o 

o 

20.  9 

36.9 

33.9 

1.77 

1.62 

25° 

28.2 

42.5 

35.8 

1.51 

1.27 

30° 

27.0 

45.4 

50.9 

1.68 

1.89 

Average  pulse  amplitude  ratios 

1.65 

1.59 

Body-wave  residuals  relative  to  0B2NV 

.21 

.20 
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Azimuth  distance  plots  of  magnitude  d i f f erent i al s 
for  the  station  RKON-NTNV 
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Azimuth  distance  plots  of  magnitude  differentials 
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Figure  43d.  Azimuth  distance  plots  of  magnitude  differentials 
for  the  station  NTNV-OB2NV 
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show  plots  of  four  categories  of  magnitude  differentials  Am.  for  various  sta- 

b 

tion  pairs  as  functions  of  azimuth  and  distance  to  NTS.  The  designations  are 
as  follows : 

A  (Am^  <  -  o  (low  apparent  attenuation  at  NTS) 

B  -  a  <  (Am^  -p)  <  0 

C  0  <  (Am^  -p)  <  +  o 

D  (Am^  >  +  o  (high  apparent  attenuation  at  NTS) 

where  p  is  the  mean  and  a  is  the  standard  deviation  of  Am^'s  for  the  given 
station  pair.  No  significant  effect  of  distance  or  azimuth  is  apparent  with 
the  exception  of  the  low-differential  "A"  events  at  1  ;  600,  6  ~  40°,  events 
from  the  Mid-Atlantic  ridge  on  the  azimuth  through  RKON.  These  events  depart 
the  source  at  a  shallower  angle  to  RKON  than  to  0B2NV  and  thus  the  RKON  wave¬ 
forms  may  suffer  more  attenuation  at  the  source  than  do  the  0B2NV  waveforms. 

Effect  of  Attenuation  on  the  m^  Measurement 

The  low  significance  of  the  "uncorrected"  magnitude  differential  between 
the  RK0N-0B2NV  station  pair  was  an  unexpected  result  of  this  study.  The  con¬ 
ventional  "corrected"  magnitude  differential  and  the  differential  in  t*  cannot 
be  distinguished  on  the  95!”  confidence  level  from  the  regional  WHS  values 
established  by  previous  studies  with  the  present  accuracy  of  measurements. 

Arguments  estimating  the  decrease  of  amplitude  of  a  single  frequency  component, 

1  Hz,  for  example,  require  a  decrease  of  amplitude  equivalent  to  .27  magnitude 
units.  The  low  significance  of  the  "uncorrected"  magnitude  differential,  how¬ 
ever,  seems  to  indicate  that  no  difference  exists  in  amplitudes  between  RKON 
and  0B2NV. 

Re-examination  of  the  method  used  in  computing  m,  reveals,  however,  that 

b 

"uncorrected"  m^,  and  conventional  m^  as  well,  is  not  a  suitable  measure  of 
the  actual  amplitudes.  The  problem  of  measurement  stems  from  instrument  cor¬ 
rection  that  assumes  the  signal  is  monochromatic,  consisting  of  a  single  fre¬ 
quency  component  with  the  dominant  period  T.  For  a  real  signal  attenuation 
changes  not  only  the  amplitude,  but  also  the  period.  The  mean  dominant  period 
of  teleseismic  arrivals  at  RKON  is  about  .7  second,  while  at  0B2NV  it  is  about 
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.9  second.  Instrument  amplification  is  higher  at  .7  second  than  at  .9  second 
by  a  factor  greater  than  1.5.  Therefore,  "correction"  for  instrument  response 
tends  to  restore  amplitude  at  0B2NV  that  has  actually  been  lost.  The  correc¬ 
tion  might  be  more  revealing  if  it  were  done  in  the  frequency  domain  using  the 
same  frequency  at  each  station.  Dividing  by  the  dominant  period  T  increases 
the  differential  for  the  conventional  since  T  is  longer  for  0B2NV. 

Observed  magnitude  differentials  are  roughly  consistent  with  the  size  of 
At*.  Given  the  subjectivity  of  picking  the  "dominant"  periods  and  amplitudes 
and  the  approximate  nature  of  "corrections",  the  noise  level  of  the  procedure 
for  computing  m^  is  similar  in  size  to  the  biases  under  study,  even  if  the 
scatter  of  body-wave  amplitude  measurements  is  disregarded.  Since  P  wave 
spectra  vary  considerably  no  exact  functional  relationships  between  the  At* 
and  the  expected  Am^  can  be  formulated.  This  reasoning  indicates  that  one 
would  be  better  off  not  using  any  period  dependent  instrument  correction. 

Even  division  by  dominant  period  is  a  questionable  practice. 

The  magnitude  differentials,  no  matter  how  poor  a  measure  of  energy  loss 

they  may  be,  can  be  causally  related  to  anelastic  attenuation,  and  are  about 

the  size  to  be  expected  based  on  At*  ~  .2.  Reciprocity  would  then  imply  that 

signals  from  NTS  would  be  weakened  by  about  the  same  amount  expressed  in  m, 

b 

units  (.18  m.u. ) . 

Travel  Time  Studies 

Although  travel  times  were  read  for  all  arrivals  that  had  reasonably 
good  S/N  ratio,  difficulties  were  encountered  during  attempts  to  find  travel¬ 
time  anomalies.  The  source  of  these  difficulties  was  inaccuracies  in  the 
locations  given  in  the  preliminary  event  lists  used.  We  concluded  that  the 
travel-time  study  had  to  be  deferred  until  better  locations  and  depths  are 
computed  using  the  maximum  number  of  available  stations.  The  ISC  list,  when 
updated,  may  provide  such  improved  locations. 

Bias  I)ue  to  Unequal  Thresholds  at  the  Various  SDCS  Stations 

Since  the  noise  levels  of  stations  used  in  this  study  vary  widely,  there 
is  a  danger  that  the  unequal  noise  levels  will  change  the  effective  bias 
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differences  measured.  The  mean  of  magnitudes  measured  near  the  noise  level 

tends  to  be  higher  in  value  than  the  actual  or  true  value  (Herrin  and  Tucker, 

1972;  von  Seggern  and  Blandford,  1976).  Although  only  events  with  good 

signal-to-noise  ratio  were  measured,  the  danger  still  exists  that  our  S/N 

criteria  will  bias  Am,  determinations.  Plots  of  Am  against  magnitude,  how- 

b  b 

ever,  do  not  show  any  clear  trend  (Figure  44),  and  we  are  reasonably  confident 
that  there  is  no  appreciable  bias  due  to  varying  thresholds.  Moreover,  our 
threshold  is  also  constrained  by  the  event  list  used,  since  we  looked  for 
already  detected  events  only. 

F.valuation  of  Attenuation  Under  the  SHOAL  Bite 

The  12.5  kt  contained  nuclear  explosion,  SHOAL,  was  detonated  in  a 
granite  body  in  western  Nevada.  The  purpose  of  the  SHOAL  project  was  to 
compare  seismic  radiation  from  this  explosion  with  earthquakes  in  the  same 
area  (Fallon  earthquakes).  The  LRSM  station  SZNV  was  operated  from  5  January 
1963  to  8  February  1963  and  was  located  about  six  hundred  feet  north  of  the 
ground  zero  of  SHOAL,  on  the  same  granite  body  where  the  nuclear  device  was 
emplaced.  The  existence  of  station  SZNV  offered  an  opportunity  to  compute 
attenuation  under  this  test  site  by  comparing  spectra  of  P  waves  observed 
here  to  those  recorded  at  a  station  on  a  stable  platform  (shield)  tvpe  struc¬ 
ture.  The  site  selected  for  comparison  is  Sleepy  Eve,  Minnesota  (SEMN),  also 
on  granite.  The  standard  procedure  described  above  was  followed  for  the 
computation  of  spectral  ratios  for  these  stations. 

Figures  El  to  E21  of  the  Appendix  show  the  vertical  SP  waveforms,  power 
spectra  (uncorrected  for  instrument  response),  and  the  amplitude  spectral 
ratios  with  straight-line  fits  for  the  21  teleseismic  events  which  had  the 
best  signal-to-noise  ratios  during  the  existence  of  SZNV.  Epicentral  data 
for  these  events  are  given  in  Table  VT.  The  signal-to-noise  ratios  of  these 
events  are  still  quite  variable  so,  due  to  the  built-in  signal-to-noise  ratio 

Herrin,  E.  and  W.  Tucker  (1972).  On  the  estimation  of  bodv-wave  magnitudes. 

Report  to  AFOSR,  Dallas  Geophysical  Observatory,  Southern  Methodist 
University,  Dallas,  Texas. 

von  Seggern,  D.  H.  and  R.  R.  Blandford  (1976a).  Seismic  threshold  determina¬ 
tion,  Bull,  Seism.  Soc .  Am.,  66,  753-788. 
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criteria  in  our  computational  procedure,  the  frequency  range  of  the  signal 
spectral  ratios  used  for  determining  At*  is  also  variable.  While  eighteen 
spectral  ratios  (SZNV/SEMN)  have  a  negative  slope,  indicating  less  high 
frequency  content  for  P  waves  at  SZNV,  only  three  spectral  ratios  have  a  posi¬ 
tive  slope.  The  relative  lack  of  high  frequencies  at  SZNV  is  quite  visible 
on  most  signals  in  the  time  domain,  where  SZNV  waveforms  have  broader  pulses 
and  less  high  frequency  riding  on  the  signals.  In  some  cases,  such  as  the 
16  January  15:09:11.4  GMT,  5  February  20:39:20.4  GMT,  and  6  February  03:27: 
56.9  GMT  events,  the  effect  is  quite  dramatic.  In  other  cases  there  is  no 
clearly  visible  difference  (16  January  05:44:54.4  GMT). 

A  notable  "counter  example"  is  the  8  January  15:46:44.1  event  that  showed 
a  broad  pulse  at  the  arrival  time  of  the  P  wave  at  SEMN.  Examination  of  film 
records  of  this  event,  however,  show  highly  unstationary,  high  amplitude, 
pulsing  long-period  noise  ahead  of  the  signal  (a  similar  pulse  is  shown  just 
ahead  of  the  signal).  Thus,  the  possibility  that  the  waveform  shown  is  a 
super-position  of  a  long-period  noise  pulse  and  the  actual  higher  frequency 
signal  cannot  be  ruled  out. 

One  possible  source  of  error  in  the  determination  of  At*  might  be  fre¬ 
quency  dependent  radiation  patterns  for  the  events  used.  However,  only  one 
of  the  events  (29  January  09:21: 16:. 9  GMT)  shows  first  motions  with  different 
polarities  at  the  two  stations.  For  the  remaining  events  both  stations  are 
on  the  same  lobe  of  the  radiation  pattern  as  indicated  by  first  motions  or  the 
general  matching  character  of  the  wave-trains.  Though  most  events  occurred 
along  the  circumpacif ic  seismic  belt,  which  can  introduce  preferred  orienta¬ 
tions  of  the  focal  mechanisms  relative  to  the  stations  used,  the  likelihood 
of  this  situation  causing  a  significant  bias  is  low  because  a  variety  of  focal 
mechanisms  exists  in  these  regions.  In  addition,  the  two  events  not  along  the 
NNW-SSE  azimuths  that  are  characteristic  of  circumpacif ic  events  show  similar 
At*  values. 

Table  VI  displays  the  At*  values  derived  from  the  spectral  ratios.  The 
mean  of  all  values  in  the  table  is  At*  =  .18  +.10  (95%  limits),  which  explains 
the  EUS  versus  WUS  magnitude  bias  of  .25  magnitude  units  at  1  Hz.  These 

i 

values  are  of  the  same  order  as  other  determinations  of  the  regional  WUS-EUS 
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magnitude  bias  and  At*  (Der  and  McElfresh,  1976;  North,  1976;  Booth,  Marshall 
and  Young,  1974). 

Observations  of  seismic  waves  from  the  nuclear  explosion  SHOAL  helps  to 
solve  the  reciprocal  problem  and  to  compute  the  attenuation  of  P  waves  to 
various  North  American  stations.  The  source  spectrum  used  for  computing 
spectral  ratios  was  obtained  by  scaling  up  the  standard  5  kt  granite  source 
spectrum  (von  Seggern  and  Blandford,  1972)  using  cube-root  scaling.  Because 
of  the  relatively  small  difference  in  yield,  the  spectrum  does  not  change 
radically,  even  if  this  scaling  law  is  not  strictly  valid  (Per  and  McElfresh, 
1976).  The  source  power  spectrum  used,  which  includes  the  LRSM  system 
response,  is  shown  at  the  top  of  Figure  46. 

Figures  45  through  47  display  seismic  traces,  power  spectra  (not  cor¬ 
rected  for  instrument  response),  and  station-to-source  amplitude  spectral 
ratios  with  least  squares  straight-line  fits.  Table  VII  provides  the  slopes 
and  the  corresponding  t*  values.  Most  of  the  values  are  above  .4;  only  one, 

HNME,  is  unexpectedly  low.  The  average  of  all  values  is  .42,  which  is  in 
fairly  good  agreement  with  other  shield-WUS  type  paths  (Per,  1976). 

The  average  .42  for  shield  to  SHOAL  site  combined  with  the  SZNV-SEMN 
differential  At*  ^  .18  yields  the  absolute  average  t*  of  .24  for  shield-to- 
shield  type  paths.  This  value  is  again  in  very  good  agreement  with  past  work 
and  other  data  from  the  literature  that  Der  (1976)  summarized. 

Spall  and  pP  can  change  this  source  spectrum.  The  measured  uphole  time 
at  SHOAL  was  about  .08  sec,  yielding  a  two-way  time  of  .16  sec  which  trans¬ 
lates  into  a  6.25  Hz  first  null  of  the  spectrum  for  total  reflection.  This 
figure  is  beyond  the  frequency  range  of  our  computations.  The  effect  of  such 

Der,  Z.  A.  and  T.  V.  McElfresh  (1976a).  Short-period  P  wave  attenuation  along 
various  paths  in  North  America  as  determined  from  P  wave  spectra  for  the 
SALMON  nuclear  explosion.  Bull.  Seism.  Soc.  Am.,  66 ,  1609-1622. 

von  Seggern,  D.  H.  and  R.  R.  Blandford  (1972).  Source-time  functions  and  spec¬ 
tra  for  underground  nuclear  explosions,  Geophys.  J.  R,  Astr.  Soc.,  31 ,  89-97, 

Der,  Z.  A.  (1976),  On  the  existence,  magnitude  and  causes  of  broad  regional 
variations  in  body-wave  amplitudes  (magnitude  bias),  SDAC-TR-76-8, 

Teledyne  Geotech,  Alexandria,  Virginia  22314. 


Figure  46.  P  wave  power  spectra  (uncorrected  for  instrument  response) 
of  P  waves  from  SHOAL.  The  estimated  source  spectrum  of 
SHOAL,  also  modified  by  the  instrument  response,  is  shown 
on  the  top  of  the  figure. 


Figure  47.  Station-to-source  amplitude  spectral  ratios  for  SH0A1 
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TABLE  VII 


Absolute  t*  Values  along  Paths  from 
The  SHOAL  Nuclear  Explosion  to  Selected  LFSM  Stations 


-log  10 

STATION  a  SLOPE  t*  =  - —  a 

7T 


NPNT 

-.657 

.48 

RKON 

548 

.40 

LVLA 

-.675 

.49 

EUAL 

-.581 

.43 

LSNH 

-.646 

.47 

HNME 

-.360 

.26 

Average 

.42 
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modification  in  spectra,  if  not  considered,  would  produce  a  slight  decrease 
of  t*  values.  Thus,  t*  values  could  be  hip, her  than  given  here. 

Vertical  acceleration  records  show  a  free-fall  time  at  the  surface 
ranging  from  .35  to  .45  seconds  at  some  of  the  close-in  stations.  Oddly, 

station  10  close  to  ground  zero  showed  little  spall.  The  free-fall  times  read 

from  the  various  accelerometer  records,  augmented  with  additional  delays 
(only  approximately  known)  stemming  from  the  two-way  travel  time  of  impulses 
between  the  shot  and  the  free  (or  spall)  surface,  yielded  delay  times  in  the 
range  of  .35  to  .65  second.  Delays  of  such  magnitude  will  introduce  at  least 
two  or  three  maxima-minima  modulation  cycles  in  the  spectra  within  the  fre¬ 
quency  range  used  here.  The  smoothing  applied  and  the  straight-line  fitting 
to  spectral  ratios  will  minimize  the  effect  of  such  modulations  on  the  attenu¬ 
ation  measurements.  In  addition,  the  sharp  fall-off  of  the  spectra  with 

frequency  cannot  be  attributed  either  to  spalling  or  pP  modulation;  thus,  it 

results  from  attenuation. 

Analysis  of  the  LRSM  Stations  EKNV  and  CUNV 

To  see  the  characteristics  of  P  waves  in  a  wider  area  of  Basin  and  Range 
near  NTS,  recordings  were  read  of  teleseismic  P  waves  at  LRSM  stations  EKNV 
and  CUNV;  RKON  was  used  for  reference  in  eastern  North  America.  While  EKNV 
was  operational  for  many  years,  CUNV  was  occupied  for  only  seven  months. 

Table  VIII  lists  the  39  events,  readings,  and  magnitudes  at  all  stations. 

Two  of  the  events  were  excluded  from  further  analysis  because  their  epicentral 
distances  are  less  than  25°.  This  section  of  the  study  is  aimed  at  establishing 
the  existence  of  magnitude  or  spectral  differences  between  RKON  and  the  two 
Basin  Range  stations.  Figure  48  and  49  show  histograms  of  the  magnitude  dif¬ 
ferentials  involved.  Both  are  significant  at  the  95%  confidence  level,  indi¬ 
cating  that  magnitudes  at  CUNV  (Am^  =  -.284  +  .187)  and  EKNV  (Am^  =  -.267  + 

.104)  are  lower  than  at  RKON.  The  "uncorrected"  magnitude  differentials  (Fig¬ 
ures  50  and  51)  are  also  statistically  significant  at  EKNV  (Am^  unC  =  -.201  + 
.116),  indicating  that  low-frequency  amplitudes  are  also  significantly  lower; 
the  RKON  dominant  periods  are  significantly  shorter  (AT  =  .133  +  .075)  than  at 
EKNV  (Figure  52).  At  CUNV  the  sign  of  the  mean  difference  in  periods  also 
indicates  longer  P  wave  periods  than  those  observed  at  RKON  (Figure  53),  but 
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TABLE  VIII 


Events  used  and  readings  at  LRSM  stations  CUNV,  EKNV,  and  RKON. 


1  15 

11  63 

21;  6:34.0 

KURILES 

STA 

CU-NV 

EK-NV 

RK-ON 

HN-8E 

DIST 
67.  K 

67.0 
70.  6 
83.8 

AN?  T 

162.5  1.10 

236. 1  0.90 

606.8  O.PO 
278.0  0.60 

LOGIO  (A/NT) 
PD 

6.97 
6.06 
6.27 

6.98 

♦  B 

LOG1  0  (A  /P )  ♦  R 
NR 

6.01 

6.87 

6.17 

6.75 

2  U 

11  63 

6:46:19.7 

CHILE 

ST  A 

cn-Rv 

EK-NV 

RK-ON 

HN-PF 

BTST 
83.9 
R4.5 
92.3 
99. 2 

$9fl  1.60 

42.9  1.50 

116.4  i.io 
°6.2  1.20 

LOGIO  («/NT) 
PR 

6.99 
6.67 
6. OP 
6.H4 

♦  P 

LOG  10  (A /N )  ♦  B 

NB 

6.18 

6.95 

6.02 

6.02 

3  11 

11  6  3 

14:39:28.9 

CALIFORNIA 

ST  A 
CU-NV 
PK-NV 

r<N-o« 

HN-PF 

fist 

9.0 

5.5 

26.7 

33.6 

A 1 P  7 

475.6  0.90 

119.9  o.70 

5 e  6 . 9  1.80 

211.6  1.40 

tor.  10  (A/PT) 

NR 

6.6  9 
4.00 
e.77 
6.77 

♦  B 

*  «» 

*  • 

O.IITTlD 

OBIT~ED 

LOGIO  (A/M)  .  9 

NP 

"  6.C5 

"  4.84 

6.62 
e  .  87 

4  19 

11  f.  3 

11:  0:69.3 

KURILES 

ST  A 
cs-nv 

SK-NV 

RK-ON 

HN-P*' 

DIST 
57.  3 

f  r 
70.  3 

M  1. 6 

AHP  ? 

126.0  1.10 

126.6  1.30 

573.  6  0.90 

729.6  0.70 

LOGIO  (A/PT) 
PS 

6.39 

6.02 

6.24 

6.06 

♦  f; 

LOG19(A/N)  ♦  4 
PR 

6.9  4 

6.14 

6.16 

5.01 

5  21 

11  63 

21:  1:35.3 

KURILES 

STA 

cn-NV 

RK-ON 

HN-NF 

DIST 

60.2 

*2.-’ 

7-.2 

ASP  T 

36.6  1.10 

113.7  0.90 

110.7  o.9C 

LOGIO  (A/«T) 
PH 

5.  17 
6.67 
6.60 

♦  P 

LOGIC  (A/«)  ♦  3 
■P 

5.21 

5.47 

6.  c  6 

6  23 

11  6  9 

16:13:11.0 

ALEUTIANS 

STA 

pK-NV 

•>K-ON 

HN-NF 

DIST 
48.6 
6’. 2 
6“.  1 

ASP  T 

40.9  0.90 

279.5  o.op 
116.P  0.30 

LOG  10  (A  /PT) 
NR 

6.04 

6.79 

6.65 

♦  P 

T.OG1  0  (A  /P )  ♦  R 
NP 

4.04 

6.74 

5.  5C 

7  3 

12  63 

23:  3:41.6 

CHILE 

STA 

C'l-NV 

FK-NV 

RK-ON 

•IN-HF 

DIS* 
79.  H 
7*.  1 
76.  3 
68.6 

1  S'* 

414.6  1.10 

23  3.7  1.10 

367.9  0.30 

32 1 • 0  1.70 

lO>!10  (A/PT) 
PB 

6.20 

o,9« 

6.09 

6.36 

♦  P 

IOG10(A/«)  ♦  q 

NP 

6.24 

6..  02 

5.  os 

6.  49 

8  10 

12  63 

16:49:47.6 

CHILE 

STA 

CII-NV 

FK-NV 

RK-ON 

01  "T 
71.F 
72.3 
7?. 4 

A  NF'  T 

109.1  0.80 

116. P  0 • 9 0 
69.6  0.8C 

T.OGIOFA/P*) 

MP 

6.62 

6.  ro 
6.72 

♦  B 

I  OGIO  1A/»1  ♦  5 

NR 

6.  4  3 
c.  66 

6.2  7 
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TABLE  VIII  (Continued) 

Events  used  and  readings  at  LRSM  stations  CUNY,  EKNV,  and  RKON. 


«  It 

12  67 

0:47:48.3 

TONGA  ISLANDS 

*  STA 

Cfl-NV 
FK-HV 
8K-OK 

DTST 
76.  U 

7fc.S 

95.4 

AXP  ? 

24.4  1.20 

131.3  1.20 

2S.0  1.00 

LOG  10  (A/9?) 
PB 

S.  14 
S.07 

S.  3t 

♦  f> 

LOG10(S/H)  ♦  B 

xb 

'.  2? 

S.es 
s.  34 

10  11 

12  S3 

17:  9:12.3 

ALEUTIANS 

ST  1 
cn-Nv 
FK-9V 
PK-ok 
HN-XF 

OTST 

as.  1 

44.4 
'P.f 
6-.  1 

;  «i  P  t 

20.4  0.70 

71.7  1.00 

8S.0  0.40 

101.*  0.90 

log  i  o  3  s  /xt) 

9P 

4.  f  r 

s.2t 

5.2' 

7.77 

♦  r 

LOG1  0  (S  /X)  ♦  p 

9p 

4.  3C 
s.  24 

4.66 

5.67 

n  19 

12  f,  ’ 

20:  .1 1:r0.  1 

ARGENTINA 

CTi 

cir-Nv 

FK-NV 

9K-ON 

!IN-«F 

nrsr 

65.  9 
'  36.4 

>1°.  a 

91.7 

t»P  T 

fip.1  1.20 
40.1  i.io 
141.0  0.OQ 
107. p  1.10 

LOG  10  (A/*-M 

■  l 

S.SQ 
S.4  1 
'.76 
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Events  used  end  readings  at  LRSM  stations  CUNY,  EKNV,  and  RKON 
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TABLE  VIII  (Continued) 

Events  used  and  readings  at  LRSM  stations  CUNV,  EKNV,  and  RKON. 
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TABLE  VIII  (Continued) 

Events  used  and  readings  at  LRSM  stations  CUNY,  EKNV,  and  RKON. 
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48.  Histogram  of  magnitude  differentials  between  stations  RKON 
and  CUNY. 


Figure  49.  Histogram  of  magnitude  differentials  between  stations  RKON 
and  EKNV. 


NUMBER 


UNCORRECTED  Amh 

Figure  50.  Histogram  of  uncorrected  magnitude  differentials  between 
stations  RKON  and  CUNV. 


Figure  51. 


Histogram  of  uncorrected  magnitude  differentials  between 
stations  RKON  and  FKNV. 
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Figure  52.  Histogram  of  P  wave  period  differentials  for  the  station  pair 
RKON-EKNV. 


Figure  53.  Histogram  of  P  wave  period  differentials  for  the  station  pair 
RKON-CUNV. 
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the  significance  of  this  result  is  lower  than  95%  with  the  sample  size  (only 
15  events)  available.  Attempts  to  determine  magnitude  differentials  at  other 
LRSM  stations  near  NTS  were  not  successful  due  either  to  the  short  operational 
life  of  these  stations  or  to  difficulties  in  reading  films  without  proper  time 
marks.  In  summary,  the  results  demonstrate  that  both  EKNV  and  CUNV  have  nega¬ 
tive  magnitude  station  effects  relative  to  RKON,  thus  conforming  to  the  general 
pattern  of  WUS  stations  and  indicating  higher  attenuation  under  WUS.  North 
(1976)  also  showed  that  EKNV  has  a  magnitude  residual  relative  to  EUS  that  is 
comparable  to  those  of  other  stations  in  the  Basin  and  Range.  Both  EKNV  and 
CUNV  are  on  well-consolidated  sediments  so  that  they  might  be  expected  to  have 
slightly  higher  amplitudes  than  a  station  on  granite  nearby. 

General  Limitations  of  the  Values  of  t* 

Much  disagreement  still  exists  in  seismological  literature  concerning 
the  proper  value  of  t*  for  short-period  waves.  While  some  investigators  use 
values  as  high  as  1.15  for  WUS-shield  paths,  others  have  pointed  out  that 
losses  along  shield-to-shield  paths  are  small  (Douglas  et  al. ,  1974),  corre¬ 
sponding  to  a  t*  of  about  .2.  Still  others  have  reported  that  attenuation  along 
certain  shield-to-shield  paths  is  not  measurable  (Filson,  personal  communica¬ 
tion),  which  corresponds  to  t*  ^  0.  These  assertions  cannot  both  be  true, 
because  taking  the  difference  of  WUS-shield  and  shield-to-shield  t*  values 
gives  a  t*  1,  resulting  in  a  huge  magnitude  differential  (about  1.5  magnitude 
units  at  1  Hz)  between  the  two  types  of  paths.  This  is  clearly  in  conflict 
with  the  observations. 

Most  discrepancies  among  these  results  are  due  to  assumptions  about  the 
source  spectra  used  to  determine  absolute  values  of  t*.  In  reality,  the  most 
unreliable  way  to  estimate  t*  is  time-domain  waveform  matching.  While  possibly 
a  good  visual  agreement  can  be  obtained  between  observed  and  synthetic  short- 
period  seismograms  using  only  the  relatively  low-frequency  components  (f  < 

1.5  Hz)  or  the  spectrum,  if  the  corner  frequency  is  below  1.5  Hz  and  ignoring 

Douglas,  A.,  J.  A.  Hudson,  P.  D.  Marshall,  and  J.  B.  Young  (1974).  Earth¬ 
quakes  that  look  like  explosions,  Geophys.  J.  R.  Astr.  Soc.,  36 ,  227-233. 
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discrepancies  at  higher  frequencies  which  determine  t*,  even  this  crude  pro¬ 
cedure  rules  out  t*  'v  1  in  most  cases  (Somerville  et  al.,  1976).  A  still 
more  suspect  procedure  is  to  match  only  the  amplitudes  of  the  first  few  cycles 
of  the  synthetics  to  the  observed  waveforms.  This  method  brings  into  play  a 
whole  set  of  unknown  factors,  such  as  the  upper  mantle  structure,  crustal 
structure,  etc.,  and  the  effect  of  t*  becomes  inseparable  from  all  of  the 
other  factors. 

If  analysis  of  t*  is  performed  in  the  frequency  domain,  then  the  range 

of  allowable  t*  is  clearly  limited  by  both  the  data  and  the  level  of  knowledge 

of  seismic  sources.  Thus,  discussions  here  are  confined  to  large  events  (m. 

b 

>  5)  because  spectra  of  such  events  are  better  characterized  in  the  short- 

period  band  (.5-4.  Hz).  Large  nuclear  explosions  have  a  corner  frequency  of 

-2 

around  1  Hz,  and  the  theoretical  spectrum  above  this  frequency  varies  as  u 
(von  Seggern  and  Blandford,  1972;  Mueller  and  Murphy,  1971).  This  estimate 
is  in  agreement  with  spectra  derived  from  reduced  displacement  potentials  used 
here  for  determining  t*  (Der  and  McElfresh,  1975).  Assuming  greater  slopes 
(Haskell,  1964)  would  lead  to  extremely  low  or  even  negative  t*  values,  thus 
contradicting  any  large  t*.  Peppin’s  (1977)  suggestion  of  an  fall-off 
rate  for  explosions  i.i  tuff  is  not  relevant  to  this  discussion  because  only 
explosions  in  hard  rock  or  salt  were  used  to  determine  t*. 

This  situation  also  exists  for  earthquakes  because  there  is  no  indication 

-2 

that  spectral  slopes  for  large  earthquakes  are  less  than  to  (Brune,  1970). 


Peppin,  W.  A.  (1977).  A  near-regional  explosion  source  model  for  tuff, 
Oeophys.  J.  R.  Astr.  Soc.,  48 ,  331-350. 

Brune,  J.  N.  (1970).  Tectonic  stress  and  the  spectra  of  seismic  shear  waves 
from  earthquakes,  J.  Geophys.  Res. ,  75 ,  4997-5009. 

Somerville,  P.  G.,  R.  A.  Wiggins,  and  R.  M.  Ellis  (1976).  Time  domain  deter¬ 
mination  of  earthquake  fault  parameters  from  short-period  P  waves,  Bull . 
Seism.  Soc.  Am,,  66,  1459-1484. 

Mueller,  R.  A.  and  J.  R.  Murphy  (1971).  Seismic  characteristics  of  under¬ 
ground  nuclear  detonations,  Bull.  Seism.  Soc.  Am.,  61,  1965. 

Haskell,  N.  A.  (1964).  Radiation  pattern  of  surface  waves  from  point  sources 
in  multi-layered  mediums,  Bull.  Seism.  Soc.  Am.,  54 ,  377-393. 
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On  the  contrary,  both  theoretical  considerations  and  observational  data  indi- 

-3 

cated  that  the  fall-off  rate  may  be  closer  to  u  (Geller,  1976;  von  Seggern 

and  Blandford,  1972).  Computing  spectral  ratios  by  utilizing  source  spectra 
-2-3 

with  either  u  or  w  fall-off  rates  will  give  t*  values  considerably  less 
than  unity  (Noponen,  1975;  von  Seggern  and  Sobel,  1976;  Sobel,  von  Seggern, 
Sweetser  and  Rivers,  1977a, b,c)  in  agreement  with  t*  determinations  from 
nuclear  explosions. 

Smaller  explosions  and  earthquakes  can  have  overall  slopes  of  less  than 

-2 

to  in  the  short-period  band,  but  in  all  cases  where  the  spectrum  has  been 
determined  (from  close-in  measurements  for  instance)  analysis  yields  t* 
values  of  the  same  magnitude  as  already  existing  studies  which  used  only  large 
high  frequency  fall-off  rates. 

Variability  of  t*  Measurements 

While  At*  measurements  appear  to  behave  consistently  for  given  types  of 
source-receiver  pairs,  it  is  important  to  have  some  idea  of  the  magnitude  of 
scatter  expected  for  such  measurements.  One  factor  causing  the  drastic  varia¬ 
tion  in  body-wave  amplitudes  is  multipathing ,  and  this  same  factor  will. 
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undoubtedly,  cause  variations  in  the  spectra  which  can  influence  At*  measure¬ 
ments.  To  determine  the  range  of  the  spectra  variation  within  a  relatively 
small  area,  the  spectra  of  P  waves  were  computed  on  a  single  sensor  of  each 
subarray  of  NORSAR.  (The  events  were  taken  from  Blandford's  [1974]  report 
because  they  were  conveniently  available  on  digital  tapes.)  The  maximum 
variation  of  P  wave  amplitudes  among  subarrays  was  about  five  to  one,  indi¬ 
cating  considerable  multipathing  and  site  effects.  If  the  area  occupied  by 
NORSAR  is  considered  a  typical  piece  of  the  Earth's  crust,  the  variability  of 
At*  can  be  assessed  by  taking  ratios  of  these  spectra.  The  reference  spectrum 
used  was  the  average  of  spectra  at  all  sensors  for  each  individual  event. 
Individual  sensor  spectra  were  divided  by  this  average  and  straight  lines 
were  fitted  to  the  ratios  to  determine  apparent  At*  values.  Since  no  real 
variations  of  t*  were  expected  across  NORSAR,  and  since  all  stations  are  on 
hard  rock  the  scatter  of  these  values  is  thought  to  be  a  result  of  multi¬ 
pathing  across  the  array.  The  procedures  used  to  compute  the  spectra  and  the 
ratios  are  identical  to  those  utilized  for  computing  our  other  t*  values. 

Epicentral  data  for  analyzed  events  are  given  in  Table  IX,  and  the  At* 

values  at  various  subarrays  are  given  in  Table  X.  A  histogram  of  all  At* 

values  is  given  in  Figure  54.  The  standard  deviation  of  these  values  is 

small,  a  .  =  .063,  giving  a  95%  limit  of  a  single  t*  measurement  of  ±.13. 
t* 

Thus,  assuming  this  scatter  typical  for  general  t*  measurements  that  are  used 
in  establishing  a  regional  difference  in  t*  for  the  shield-to-shield  and 
shield-to-WUS  paths  only  a  few  pairs  of  measurements  should  suffice,  but 
individual  measurements  can  still  considerably  overlap.  On  the  other  hand, 
a  much  larger  number  of  amplitude  measurements  is  needed  to  establish  a  magni¬ 
tude  station  bias  for  a  pair  of  stations.  These  facts  are  quite  obvious 
from  the  above  results  for  both  At*  and  Am  . 


Blandford,  R.  R.  (1974).  Short-period  signal-to-noise  ratio  at  NORSAR, 
SDAC-TR-74-13 ,  Teledyne  Geotech,  Alexandria,  Virginia  22314. 
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TABLE  IX 
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Summary  of  At*  Values  Computed  at  NORSAR 
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N  =  139 
MEAN  =  .024 
STD.  DEV.  =  .063 


At*  values  between  various  subarrays  of  NORSAR 


DISCUSSION  OF  RESULTS 

The  purpose  of  this  study  is  the  evaluation  of  attenuation  under  NTS. 
Attenuation  is  the  only  truly  reciprocal  factor  with  respect  to  source- 
receiver  postions.  There  is  an  additional  problem  for  a  nuclear  explosion 
source:  reciprocity  does  not  hold  because  the  equivalent  elastic  source 

changes  with  shot  medium,  even  if  reciprocity  would  be  strictly  valid  for  the 
same  equivalent  elastic  sources. 

Evaluating  attenuation  characteristics  using  teleseismic  P  waves  is  a 
difficult  task,  utilizing  several  basic  assumptions  that  cannot  alwavs  be 
easily  substantiated.  The  difficulty  lies  in  deciding,  which  observed 
anomalies  are  reciprocal  with  respect  to  interchange  of  sources  and  receivers. 
Ideally,  a  station  located  over  an  attenuating  upper  mantle  should  show  lower 
amplitudes  for  teleseismic  waves  than  a  station  with  less  attenuation  under¬ 
neath,  and  their  spectral  ratios  should  show  the  degree  of  differential 
attenuation.  Crustal  structure,  however,  can  change  the  absolute  level  of 
P  wave  amplitudes,  especially  if  the  surface  is  covered  with  low-velocity 
materials.  While  crustal  structure  can  also  change  the  slopes  of  spectral 
ratios,  our  calculations  using  a  large  number  of  crustal  LRSM  station  responses 
show  that  in  most  cases  the  slopes  do  not  change  significantly  in  the  total 
frequency  band  (.5-5.  Hz)  even  if  the  absolute  levels  of  P  wave  amplitudes 
are  modified  significantly  by  the  receiver  structure.  A  similar  argument  can 
be  used  for  the  effects  of  systematic  crustal  focusing  or  defocusing:  multi- 
pathing  could  probably  be  compensated  for  by  using  many  events  at  various 
azimuths  and  distances.  In  addition,  the  possibility  exists  that  at  some 
individual  stations  the  crustal  response  might  significantly  change  the  shape 
of  observed  body-wave  spectra  so  that  t*  cannot  be  reliably  measured.  Still, 
the  spectral  ratio  measurements  remain  more  reliable  diagnostics  of  attenuation 
and  thus  of  the  expected  effect  on  m^  from  sources  at  the  station  site  than 
arc  absolute  amplitudes. 

Measuring  attenuation  under  individual  stations  is  more  difficult  than 
>■  ;!ng  regional  attenuation  patterns  because  (Per,  1976;  Per  and  McElfresh, 

•  •' )  regional  studies  average  over  many  stations  and  tend  to  minimize  the 
i.  r  cal  crustal  distortion  of  amplitudes  and  spectra  at  individual 
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An  apparent  discrepancy  exists  between  the  sizes  of  magnitude  and  ampli¬ 
tude  differentials  at  RKON  and  0B2NV  on  one  hand  and,  on  the  other  hand,  the 
size  of  At*  .17.  This  is  possibly  due  to  crustal  focussing  at  NTS.  The 
magnitude  differential  is  significantly  less  than  predicted  by  the  amplitude 
attenuation  factor  computed  for  a  single  frequency  (1  Hz  was  commonly  used  in 
previous  studies). 

The  magnitudes  are  also  significantly  lower  by  a  similar  amount  at  NTNV 
and  NT2NV  relative  to  RKON,  if  the  .2  magnitude  unit  correction  for  parallel 
layer  crustal  effects  mentioned  above  is  made.  This  correction  assumes  that 
the  crust  and  the  upper  mantle  under  the  volcanics  at  Pahute  Mesa  and  Climax 
Stock  are  the  same,  which  may  not  be  true. 

Based  upon  measurements  of  dominant  P  wave  period  differentials  between 
EUS  and  WUS  (.1-.2  sec),  it  is  unlikely  that  the  previously  reported  EUS-WUS 
magnitude  differentials  (Booth,  Marshall  and  Young,  1974)  can  be  attributed 
solely  to  changes  in  period.  Even  if  a  .1  magnitude  unit  (maximum)  difference 
is  assumed  as  result  of  the  period  differential,  the  remainder  must  be  due  to 
differences  in  absolute  amplitude. 

The  study  of  SHOAL  and  the  comparison  of  SZNV  and  SEMN  indicates  that  the 
mantle  under  SHOAL  is  typical  of  the  WUS  both  in  the  EUS-WUS  differential  At* 
and  in  the  absolute  values  of  t*  for  SHOAL-to-shield  type  paths.  The  magnitude 
bias  at  the  SHOAL  site  has  been  shown  to  be  not  significantly  different  from 
MNNV  (Der  and  McElfresh,  1976),  which  had  a  bias  of  -.28  +.05  relative  to  glo¬ 
bal  averages  of  m^  (Booth,  Marshall  and  Young,  1974). 

The  study  of  CUNV  and  EKNV  shows  that  these  stations  are  also  typical  of 
WUS  with  significantly  lower  magnitudes  and  larger  P  wave  periods  than  EUS 
stations. 

This  study  firmly  establishes  that  short-period  P  waves,  crossing  the 
upper  mantle  under  the  Basin  and  Range  province,  lose  a  great  portion  of  their 
high-frequency  energy  relative  to  those  travelling  under  shields.  The  number 
of  spectral  measurements  presented  here  is  greater  than  in  any  similar  data 
set  previously  published.  The  differences  in  spectral  content  are  statisti¬ 
cally  significant  on  at  least  the  95%  confidence  level. 
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There  seems  to  be  the  possibility  that  at  individual  stations  the  measured 
magnitude  bias  of  incoming  teleseismic  waves  is  significantly  modified  at  NTS 
by  crustal  structure  resulting  in  focussing.  Such  modification  is  not  in  gen¬ 
eral  reciprocal  with  respect  to  exchange  of  source-receiver  positions,  while 
the  effect  of  attenuation  is  reciprocal. 

The  spectral  ratio  results  presented  in  this  report  fit  well  into  the 
framework  of  previous  studies  of  t*  in  shield  regions  versus  shield-WUS  type 
paths.  Figure  55  shows  a  compilation  of  t*  calculations  plotted  against  epi- 
central  distance  previously  presented  in  Seismic  Data  Analysis  Center  Techni¬ 
cal  Report  76-8  with  some  new  data  points  added.  On  the  right  of  the  figure 
we  show  the  EUS-WUS  coverage  differential  in  t*  for  some  Novaya  Zemlya  shots 
and  for  the  two  station  pairs  SEMN-SZNV  and  0B2NV-RK0N  presented  in  this 
report.  The  differentials  in  t*  are  roughly  equal  to  differences  between  the 
absolute  t*  determinations  at  teleseismic  distances.  Thus,  within  the  limits 
of  accuracy  of  the  present  measurements,  the  NTS  stations  behave  as  typical 
WUS  stations  with  respect  to  t*. 
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Figure  55.  Summary  of  available  t*  calculations  vs.  epicentral  distance. 
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